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Abstract 

We generalize the earlier theory by Predrickson [J. Rheol. 38, 1045 (1994)] to 
study the orientational behaviour of the hexagonal phase of diblock copoly- 
mer melt subjected to steady shear flow. We use symmetry arguments to 
show that the orientational ordering in the hexagonal phase is a much weaker 
effect than in the lamellae. We predict the parallel orientation to be stable 
at low and the perpendicular orientation at high shear rates. Our analy- 
sis reproduces the experimental results by Tepe et al. [Macromolecules 28, 
3008 (1995)] and explains the difficulties in experimental observation of the 
different orientations in the hexagonal phase. 

I. INTRODUCTION 

Polymeric liquids subjected to shear flow demonstrate a very peculiar phase behaviour. 
Their phase diagrams not only contain regions of stability of the different symmetry types 
(lamellae, hexagonally packed cylinders and so on), but also these regions have an inter- 
nal structure. Application of shear breaks the rotational symmetry selecting the preferable 
direction. Thus, at given external parameters (temperature and shear rate) a certain ori- 
entation of the symmetry pattern with respect to the selected direction is more stable than 
the others IfO!. 
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Experimentally this phenomena was first observed for the lamellar phase 0. In this 
work a microphase separated poly(ethylene-propylene)-poly(ethylethylene) diblock copoly- 
mer melt was subjected to an oscillatory shear. Near the order-disorder transition (ODT) 
the lamellae have their normal parallel to the velocity gradient (the parallel orientation) at 
low shear frequencies, while at high shear frequencies the lamellae have their normal par- 
allel to the vorticity direction of the shear flow (the perpendicular orientation). At lower 
temperatures the parallel orientation is always the most stable one. 

The first theoretical attempt to study the orientational phase transition was by Cates 
and Milner ||. They considered the equation of motion for the order parameter with a 
coupling of applied steady shear to the composition fluctuations. They found that flow 
completely changes the fluctuation spectrum. As a result, the fluctuations of the order 
parameter are suppressed and the ODT temperature is raised. The authors suggested that 
the perpendicular orientation is stable near the ODT since the composition fluctuations with 
wave vectors normal to both the velocity and the velocity gradient are the least affected by 
shear. 

Later Fredrickson has shown that the angle-dependence of the fourth order vertex func- 
tion is crucial for construction of a realistic non-equilibrium phase diagram. Within this 
framework Fredrickson reproduced the experimental observation by Koppi et al. ||. The 
theory fails to describe the selection of the parallel orientation at higher frequencies Q. 

In the present paper we generalize the Fredrickson theory to describe the orientational 
phase transitions in the hexagonal phase of a diblock copolymer melt subjected to steady 
simple shear flow. The symmetry pattern of this phase consists of hexagonally packed 
cylinders made up of blocks of one type immersed in the surroundings of the other blocks (see 
Fig.|l]). In equilibrium this phase appears in between lamellae and body-centered cubic phase 
IIU. The experiments under shear flow show an analogy with the orientational behaviour 
of the lamellar phase fTJHE]- m s l° w fl° w the symmetry pattern has "2 dots up", while in the 
faster flow the "1 dot up" orientation appears (see Fig.|2|). In order to keep universality we call 
them the parallel and the perpendicular orientations, respectively. There have been many 
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observations of the ordering of the hexagonal phase subjected to shear flow 0-0. Recently, 
Tepe et al. have varied temperature and shear rate and observed both orientations of 
the hexagonal pattern in a non-symmetric polyethylene-poly(ethylenepropylene) (PE-PEP) 
diblock copolymer melt. We schematically present their results in Fig.[| Our goal is to 
reproduce this dynamical phase diagram theoretically. 

We start with a speculative analogy with lamellae. Both orientations of the hexagonal 
pattern could be considered as parallel lamellae with different interlamellar distances (see 
Fig||). High shear squeezes the pattern and the lamellae with the smallest interlayer distance 
are favorable. In very slow shear the intermolecular forces resulting in microphase separation 
will play the predominant role and lamellae with the biggest interlayer distance are stable. 
Such simple reasoning already reproduces the main features of Fig.^|. 

We proceed with more exact and motivated analysis. In Section II we formulate a dy- 
namical model and study its equilibrium limit. In Sections III and IV we apply the methods 
of Fredrickson [I]] to analyze the high-shear and low-shear behaviour. In the Conclusions we 
compare the orientational behaviour of the hexagonal and lamellar phases and explain the 
origin of difficulties in experimental study of this behaviour. 

II. DYNAMIC EQUATIONS 

Let us consider a diblock copolymer melt. The starting point of our analysis is the 
dynamic equation for an order parameter ip(r), which we choose to be a deviation of a local 
density of monomers of one block from its average value || . In the present work we use the 
Fokker-Planck equation for an incompressible block copolymer melt |I|j3]JT0[]: 

9P ri n f S 

where H[ip] is the Landau-Ginzburg Hamiltonian 
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and n is the Onsager mobility coefficient, which we assume to be constant (see [1,11,12] for 
discussion) . 

For a certain configuration of the field ip(q), the function P[ip] gives the probability 
of its realizations. In the stationary state, P[tp] reduces to the Boltzmann distribution 
P ~ exp (—H[$\) in the limit D = |T3 . 



In the derivation of the equation (|2.1| ) we assumed the standard flow geometry v = D y e x . 



It means that we ignored any alteration of the velocity profile caused by the internal structure 
of the melt. Fredrickson has showed Q that by taking into account different viscosities of 
the two blocks one can approximate the real velocity profile by the same functional form 
v = D c ffye x . Here D e s is a renormalized shear rate: 

D c s = D(l + corrections) (2.3) 

In lamellar phase at high shear, the lamellae are perpendicular at high temperatures, 
and parallel at low temperatures. Fredrickson has shown that the orientation-dependent 
hydrodynamic corrections are needed to explain this transition. In hexagonal phase such 
an effect at high shear is not known, and we postpone hydrodynamic corrections to future 
analysis. Here we approximate D e g by D. 

We take the first and the second cumulants to obtain the equations for 

C(p) = Mp)> (2.4) 

%) = (^(pW-p)}-Wp))(^(-p)) 

The resulting equations can be greatly simplified if we use the principle-harmonic approxi- 
mation for the mean density profile: 

n 

C(p) = ^ai[^,,„„« + (5 Pr9on (i)] (2.5) 
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where the set of vectors nW determine the lattice symmetry. The equation (|2.5|) says that 
all structures in the system have the same typical size ~ 1/go, where go corresponds to the 
primary peak of the structure factor S*(p) (see the first term in ( |2.2|) ). 
For the steady state the equations transform to: 

hi = rat + %+ l -a { J A(g n«, -g n«, q, -q)S(q) + ^ (2.6) 
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where we defined the structure constants: 

** = \f I e(-9bn (0 ,qi,cb)C?(qi)C7(q 3 ) 
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C =^A(p,q,-p,-q)C(q)C(-q) 

In equation ( |2.q* ) we introduced artificial external fields hi. This is equivalent to an 
introduction of the additional term in the Hamiltonian H ext = — J h(q)ip(— q), h(q) = 
Y^i=i M^qgonC*) + ^q,- go n( 4 )]) which describes the interaction with an external field h(q). The 
fields hi will allow us to construct a potential $ governing the dynamics. In equilibrium 
the potential <£> has the meaning of the free energy of the system. Thus, by introducing the 
fields hi we obtain an analytic continuation of the free energy to the dynamic case ffl. This 



is possible because our model is conservative [14 



As our last simplification we introduce an approximation for the vertex functions £ and A. 
Since we use the principle-harmonic approximation in ( |2.5| ), we assume that all wavevectors 
have the same modulus |q| = go- Moreover, following p5|-4T^1, we take into account the weak 
angle-dependence in the 4th-order vertex function A, so 

C(qi, <12, qs) = f <Kqi + q2 + qs) (2.9) 
A(qi, q 2 , -qi, -q 2 ) = A [l - (3 ( Ql ■ q 2 ) 2 ] (2.10) 

/?< 1 



where q = q/g denotes the unit vector in the direction of q. 

To make the general equations ( |2.6| ) and ( |2.7| ) specific for the hexagonal phase one needs 
to calculate the structure constants ( |2.8| ) taking into account the symmetry of the phase. 
The average density profile for the hexagonally packed cylinders is given by ( |2.5| ) , where the 
basis vectors are: 



= {0, cos0, sin0} 

= -{0, — cos0 — v^3 sin0, v^cos^ — sin0} 
= -{0, — cos0 + V^sin^, — cos0 — sin0} 



(2.11) 



The angle defines the orientation of the hexagonal pattern. The case = corresponds 
to the parallel orientation in real space (see Fig.Q). 

Using (|2.8|), ( |2.11| ) and (|2.9|), ( |2.10| ) we obtain for the structure constants: 
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% = ^a k ai 

£ = A [al (1 - 0(p • n«) 2 ) + a 2 (l - (3(p • n^) 2 ) + a 2 (l - 0(p • n®) 2 )] 

z ^ jfc ^ Z 
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We introduce the notation: 



A 



*(P) = \ J ^(q) [1 - /3(P • q) 2 ] 
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p.e.p = r + a(p) + A^a 2 [l - /3(p • n^) 2 ] 



i=i 



•Vp) = r - P • e • P+ (P- 5o) 
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and rewrite the steady-state equations in the final form: 

fc. = ( r _ n« • e ■ n«)a, + £a fc a, - ^Aa?(l - /3) 



(2.13) 

(2.14) 
(2.15) 



(2.16) 
(2.17) 
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Notation ( |2.13| )- fl2.15|) has a clear physical meaning. The fluctuation integral ( [2.13|) takes 



into account the fluctuations of the order parameter [JTHJ] and renormalizes the temperature 
in the system. Because of the angle-dependence of the 4th-order vertex function A (eq. [2.10| ), 
the renormalized temperature also has an angle-dependence. Expanding it to the first order 
in (3 (one can easily check that the anisotropy tensor is of order 0(f3)), we extract this 



angle-dependence and get (|2.14j ), where r denotes the p-independent part of the renormalized 
temperature and — p - e p adsorbs the other terms. Finally, 5o(p) is an equilibrium structure 
factor, which in the limit (3 = reduces to the one studied by Brazovskii ||18|| , Fredrickson 
and Helfand [HI. 



Equilibrium analysis 



The equations fl2.16|) and fl2.17| ) describe the behaviour of the hexagonal phase at any 
shear rate. For the particular case (D — 0, /3 — 0) these equations were studied in a number 
of articles P JT^ , [18| . Before studying the dynamics we want to show the influence of the 
angle-dependence in A (|2.10| ) on the equilibrium phase behaviour, so we construct the free 
energy at equilibrium D = 0. The equation (|2.16|) reads: 



S(p) = S (p) = r-p e -p + ip-qo) 



Performing the integration in ( p.!3| ) up to the first order in /3 we get: 



(2.18) 
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Here en denotes the trace of the anisotropy tensor e. Separating p-dependent terms in cr(p) 
with the help of (|2.14 ) we obtain: 



aX aX 2 2 

t + — + -^/5 e a + K a i + a 2 + a 



OiX (3 2 m 2 (2 ) (2) 2 (3) (3)i 

fkj = 3 Oy + PKa>i n l n) + a 2 n\ n) + a 3 n\ n) >\ 



(2.20) 
(2.21) 



The equations of motion for aj ( 2.17|) have a potential form: 



To integrate Q2.22Q one needs to take into account dependence of r on and treat an integral 



of a composite function: 



doi{...)= I dr^(...) (2.23) 



JrO 



where the Jacobian can be found from ( 2.20|) , and r$ is a temperature at which the 



disordered phase loses stability, given by equation fl2.20|) with a\ = a 2 = a 3 = 0. Integrating 
( P- 22] ) we obtain the free energy: 

$= & + $ sir (2.24) 
*' = + «(1 - § )(Vr - y/rd + ^t(3(^= - -L) + i« 2 A/?(i - 1) (2.25) 



1 1 
$ str = --A(l + p)(cq + a\ + af) + 2^aia 2 a 3 - ~/3A (a 2 ^ + a 2 a 2 + a\a%) (2.26) 

where we have separated terms determined by the structure constants. 

The amplitudes that minimize the potential $ are found by solving the set of equations 
hi = 0. Its most important feature is that a does not depend on the n®. Thus, the resulting 
equations are isotropic with respect to and the solution has the form a± = a 2 = a 3 = a. 

To explore the stability region of the hexagonal phase we need to solve the resulting 
system of equations for given temperature and composition: 

aX (aX) 2 . aX 2 2n „. 2 (n 
r = r + — + + — -a 2 /? + 3Aa 2 (2.27) 

•Jt or zr I 



aX (aA) 2 



r = r + — + ^-f /3 (2.28) 
r-i^-Aa 2 (^ + /3)+ea = (2.29) 
$= $' - ^Aa 4 (l + 2/3) + 2£a 3 (2.30) 

The calculation technique of Fredrickson and Helfand || allows us to demonstrate the in- 
fluence of the /3-terms on the equilibrium phase diagram. In Fig.^ and Fig.|5] we show that 
the approximation ( |2.1C| ) results in a non-significant quantitative shift of values. The effect 
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of (5 on the equilibrium phase diagram is tiny, but in the next section we will demonstrate 
that f3 is dominant in explaining the dynamical orientational ordering. 



III. STRONG-SHEAR BEHAVIOUR 



To study the orientational dynamics in the strong-shear regime we need to solve ( [2.16| ) 
and C P . 1 71) in the limit D — > oo. In this case it is impossible to obtain a solution of Q2.16|) 
as a perturbation series in l/D. Instead, we use approach developed by Cates and Milner 
@. One can apply the RG methods [10] to find the asymptotic behaviour of the structure 
factor S(p,D — > oo) = S^p). Then, at large shear rates the structure factor can be 
approximated by interpolating between So(p) and Soo(p). Cates and Milner introduced the 
following approximation: 

fi^p) = S \ P ) + ^(p) = r- p.e.p + (p- g ) 2 + 1 (^f) ^ (3.1) 

where Co = |(487r) 1//3 r (|). To the leading order in D the fluctuation integral is equal to [|TJ: 



a(p) = {aXf^ih - (3{I 2 pl + hp 2 y + hf z )\ 



(3.2) 



where 



and 



(48vr) 1 /3 ^ d 
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D* = /iAya 



27rr(I) 27rr(f) 2nT (f ) 



For the external fields hi the equation ( |2.17|) results in: 



- + a(n" , )+ AA(l-^/3 



o< + ^aaO^ - ^Aa 3 ( 1 + -/? 



where 
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(3.3) 



/t 2,2,2 

A = a x + a 2 + a 3 
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The potential $ is given by ( |2.22| ). Unlike the equilibrium situation we do not have any 



problem with integration ((|3.3|) does not depend on both r and e^), and the equation for 
the free energy is straightforward: 

3 -t / i \ / \ 3 



a 4 (3.4) 



$ = tA + J2 ^)a) + \\ (l - \p) A 2 + 2^ ai a 2 a 3 - ~A (l + ^ ]T , 

3=1 ^ ' ^ ' J=l 

The amplitudes aj are the solutions of the equations 

/ii = (3.5) 

In the absence of the term these equations have a uniform solution a! = a 2 = 

a 3 . The presence of this term breaks the symmetry of the equations and allows for the 
nonequal amplitudes. We construct the solutions of ( |3.5| ) as a perturbation series in two 
small parameters: D and (3. The highest order in the perturbation series for the amplitudes 
is determined by the following argument. Our goal is to obtain an angle- dependent free 
energy. This angle-dependence can only appear via various but symmetric combinations of 
the basis vectors n«. Let us discuss an example: 

n (l) 2 , +n (2)2, +n (3)2, 

This combination is angle-dependent if a; > 3. One can easily check that it is also true for the 
other symmetric combinations. The only way for such a combination to enter the equation 
for the free energy is via a term like a 3 , which is proportional to /3 3 ® 3 . Thus, in what follows 
we keep terms up to 0(/3 3 3) 3 ). Solving (|3.3|) for the amplitudes and substituting them into 
(B.4I), we obtain 



$ = $ + + $ 2 S) 2 + $ 3 D 3 (3.6) 

The coefficients $; are given by: 

-2 + 15(2 - 5x)x - 2(1 - 10xf/ 2 £ 4 
~ 125 A 3 

2 e 2 « 2/3 



$! = —h (1 + VI - 10x) 
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$ 2 = + - , 

5 11 VI - 10x 
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J? 



3 3 VT£Tto(122 - 5x) + 27(14 - 5x) 



COS( 



T 



(1 - 10x) 3 / 2 v J 7 ^ 16(8 + x) 3 

Here we have introduced the dimensionless temperature x = rA/£ 2 and kept only the leading 
terms in f3 for a given order in 35. The second term in $3 is the angle-dependent contribution 
we were looking for. 

The spinodal temperature is obtained from the condition: $ \ T=Ts = 0. It reads: 



rs = r (°) +r WX) + r f £2 (3)£3 



(3.7) 



where: 
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(«A) 2 / 3 

9/? 2 (5-2/3) 2 (I 3 -I 2 ) 2 a 4 / 3 A 7 / 3 
64(13 - 7/3) £ 2 
81(/ 3 -/ 2 ) 3 /3 3 (5-2/3) 5 /«A 2 ^ 2 



COS I 



512(13-7/3) 3 V. £ 2 
Finally, we calculate the transition temperature between cylinders and lamellae. Up to 



the first order in D we obtain: 



rtr = -I±M^ + J l(7 + 3V6)(a\f 3 ® 
5 A 



(3.8) 



Now we summarize our results. In the high-shear limit the hexagonal phase is found 
to be stable in the temperature range from t s to r tr (see Fig.|]). Here we suppose that 
hexagonal phase is the first phase appearing at cooling down from melt. If this is not a case, 
the temperature range (r s ,r tr ) transforms to (r*,T tr ), where is a transition temperature 

from a hypothetic (bcc, gyroid, ...) to the hexagonal phase. The factor in front of cos 60 

(3) 

in the equation for r s is negative. This means that the spinodal temperature is higher 



for the orientation with </> = n/6. Moreover, the free-energy is minimal for this orientation 
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for all values of temperature from the range (r s ,r tr ). Thus we predict the perpendicular 
orientation to be the only stable orientation in the high-shear limit. 

We emphasize that the appearance of the angle-dependence in the free energy only in 
the 0(/3 3 D 3 ) is not a coincidence. It reflects the internal symmetry of the system. Any 
figure on plane can be oriented with respect to the particular direction only if its shape 
deviates from circle. In other words, if one expands the figure's shape into the plane waves 
around circle: r(0) = Rq[1 + ^2 a n e m ^}, the interaction with a selected direction will appear 
in the first non-zero order. For hexagons it gives n = 3. The n-th order in the expansion 
corresponds to the interaction between n different wavevectors. Taking into account that the 
parameter (3 introduces the pare-wise interaction between vectors (see ( |2.10| )), we conclude 
that the angle-dependent terms in the free energy should be at least 0(/3 3 ). The role of the 
parameter 2) is different. Shear breaks the rotational symmetry in the system, stipulating 
the preferred (gradient) direction and allowing for a discrimination between orientations. 
Thus the angle-dependent terms in the free energy should be proportional at least to the 
lowest possible power of D. This lowest possible power is determined by the concrete form 
of the structure factor S(k). In ( |3.2| ) (3 and D enter as a one combination. So, one would 
expect the free energy to be angle-dependent starting from the 0(f3 3 ® 3 ) order. 

We illustrate the concept with the following example in which we calculate the free energy 
of systems with other rotational symmetry. We consider for the moment an artificial phase 
which the average density is given by: 

2 

C( k ) = Yl a *( 5 k, 90 n« + <*k ,-gonW) (3-9) 
i=l 

where rij are two vectors lying in the plane perpendicular to the velocity direction, -n^ = 
cos#. We consider two ) 9 = 7r/2, which describes square packed cylinders and 

b) 9 = 7r/12, which describes highly non-symmetric cylindrical pattern. Calculating the 
structure constants (|2.12 ) and integrating (|2.17|) with the help of ( |3.2|) , we obtain for the 
free energies: 
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4(4 + v / 3)/i + 12(/ 2 + / 3 ) 
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;2 + v / 3) 



cos . 
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(3.11) 



In the case a) the structure deviates from a circle in the second order (a = a.\ = , «2 7^ 0), 
while in the case b) already in the first order. This gives (3 2 and (3 standing in ( |3.10|) and 
( [3. 11| ) in front of the angle- dependent terms. Both cases have the same power of 53 as (3 
because of ( |3.2|) . In the next section we shall see how a different expression for the fluctuation 
integral will produce the different lowest possible power of D. 



IV. WEAK-SHEAR BEHAVIOUR 



In this section we consider the other limit D — » 0. In this case the solution of ( |2.16| )- 
( |2.17| ) only slightly deviates from the equilibrium one and we can construct a perturbation 
theory with a small parameter D. Thus the structure factor is given by the equation: 

d-i n 



00 , v n 



PxS (p) 



d P y 



So(p) 



(4.1) 



where S (p) is the equilibrium structure factor. 

Following the conclusions of the previous section we evaluate S(p) up to 0(/3 3 ) and the 
lowest possible order in D, which is 0(D). However, S'^^p) does not contribute to the 
fluctuation integral. Therefore, we keep S^ 2 \p) in the expression for the structure factor. 
Performing integration in ( |2.13| ) with 5*(p) given by ( |4.1| ) and using ( |2.15| ) for S*o(p) we 
obtain the fluctuation integral: 
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and cr eg is an extension of (|2.19| ) to the higher orders in (5. Here we assumed the summation 
over repeated indices. 

In the equation ( |4.2|) and in what follows the terms of 0(/? 3 ) are cumbersome and we do 
not present them here. However, we did use them in our calculations. 
The p-independent terms in (4.2) contribute to the equation for r: 



aX 



r = r + XA + - r R 1 -l- 



D V (aX)'' 
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For the anisotropy tensor (|4.2|) together with ( (2.14T ) gives a closed set of equations. 
Iterating them up to 0(D 2 ) and 0(/3 3 ) we obtain: 

£> V (c*A) 3 7r p(a) 



e • = - I — I v "' v " E lZj 4- 
CiJ *J \ d* ) 24r 7 /2 « 



(4.4) 



where 
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2 '^«+3^E^h+2«rx* , (2-%) 

s=l 



+ ■ ■ 



E-f = j {35 ix 5 jx + 35 iy 5 jy + S iz 5 jz } + f3 2 j ^^§75 [59(4^ + 5^5^) + 225 iz <^] + 
A 3 

^ + 2n« 2 (<f fa <f i!B + 25 iy 5 jy ) + 4n? ) 7i$' ) (l - <%) 
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Equations ( |4.3|) and ( [4.4|) are enough to construct the free energy $ as a solution of 
( |2.22| ), where hi, r and are given by ( |2.17| ), ( f4.3| ) and ( |4.4j ), respectively. Integration of 
( |2.22| ) using ( |2.23| ) leads to the following equation for the free energy: 

Z = F -fi(£-\\< 3 -?^ ms64> (4.5) 



\D*J 1280r 

where a is the equilibrium amplitude given by the equations (|2.29|) and ( p.27| ). In (|4.5|) F 
stands for the angle-independent part of the free energy. 
The spinodal temperature r s is given by: 

^ t + i33 (^) 2 ° 4 S C086 * (4 ' 6) 

where t includes the angle-independent terms. The spinodal temperature is maximal for 
= 0. At the same time, the free energy is minimal for the same orientation. It means that 
only the parallel orientation is stable under low shear. 

To finish our symmetry analysis we note that in principle the fluctuation integral ( [I.2| ) 
contains all possible combinations of powers of (3 and D 2 . According to the symmetry 
arguments the angle-dependence appears in the lowest possible order in D (which is D 2 in 
this case) and the third order in f3. This is in agreement with fl4.5|) and G4.6|) . 



V. CONCLUSIONS 

In the present paper we have shown how the orientational behaviour of the hexagonal 
phase under simple shear flow can be described in the framework of the dynamical model 
first developed by Fredrickson for the lamellae 0. In this model the angle- dependence of the 
4th order vertex function A (see P-10| ) plays a crucial role, although it is of no importance in 
equilibrium (see Section II for details). The parameter (3 introduces the interaction between 
structure and shear flow and allows for the rotational symmetry breaking. The character of 
the interaction depends on the shear rate and the resulting phase diagram has a complex 
structure. We predict the parallel orientation to be stable at low and the perpendicular ori- 
entation at high shear rates. Our results are in agreement with the experimental dynamical 
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phase diagram (FigfJ) for the PE-PEP system |7j. Our analysis shows that experimental 
study of the hexagonal phase orientation is actually very difficult to perform. Because of 
the symmetry arguments, the difference in the free energies of the different orientations 
A$ = $x — $11 is proportional to 



^3D! ) D ^ ^ 



This difference is extremely small (compare with A$ oc (3 for lamellae and a particu- 
lar orientation is only slightly more stable than the other one. Some crude experimental 
by-effects (like the ordering influence of the sample walls) can easily suppress the pure ori- 
entational behaviour we discussed. 

In comparison with the work [|TJ we ignored the hydrodynamic corrections to the effective 
shear rate _D e ff (see Eq.(2.27) in [Q). These corrections are responsible for the J_— >|| transition 
in the lamellar phase under strong shear. Taken into account they would probably lead to 
the same transition in the hexagonal phase. However, this is not an easy task because one 
should keep terms up to 0(/3 3 ® 3 ) in order to have an angle- dependent value. 

Another remaining problem concerns the stability of transversal cylinders, i.e. cylinders 
with their axis oriented not in the direction of flow [|]]. Our model cannot describe stability 
of the transversal orientation ||. Another approach should be developed to complete our 
study. 
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FIGURES 




FIG. 1. The molecular structure of a cylinder 



FIG. 2. Schematic dynamical phase diagram of hexagonal phase under shear flow: the parallel 



0) and perpendicular (0 = 7r/6) orientations (see eq.2.11 
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FIG. 4. The spinodal temperature (solid line: (3 = 0, dashed line: [3 = 0.5) 
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FIG. 5. The transition temperature between lamellae and cylinders (solid line: (3 = 0, dashed 
line: (3 = 0.5) 




FIG. 6. Region of stability of the hexagonal phase in high shear flow 
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